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Recently Gasser et al. [QQ have presented updated 
conversion factors between gamma radiation air kerma 
and concentrations of 40 K and the 226 Ra and 232 Th se¬ 
ries radionuclides within soil. Such conversion factors 
are important for calculating external exposures due to 
naturally occurring radioactive elements in the earth’s 
crust [2- The authors derived the conversion factors 
using MCNPX Monte Carlo simulations. The reported 
conversion factors are on average 20% smaller than es¬ 
tablished estimates from both experiment )3J and simu¬ 
lation [4] (Table 5 of their paper). They attributed the 
difference as mainly due to their use of more up to date 
emission probabilities and branching ratios for radioac¬ 
tive decay. 

In our view the predominant cause of the difference is 
in fact a finite size issue with their simulations. Their 
simulation setup consisted of a cylindrical source con¬ 
taining radionuclides homogeneously within soil down to 
a depth of 1 m, and with varying radius up to 35 m. The 
air kerma was calculated for 1 m above the ground sur¬ 
face. From Fig. 3 in their paper they concluded that a 
15 m radius source cylinder is sufficiently large for the 
conversion factors to reach their asymptotic values for a 
true half-space geometry. The infinite half-space is the 
applicable geometry for the conversion factors defined by 
convention in ICRU publication 53 [5]. 

A 35 m radius simulation is not sufficient to model a 
half-space, as can be seen from Fig. [l] in this letter. The 
figure displays some of our results covering a much larger 
range of source radii than considered by Gasser et al. 
Our results were calculated using the PHITS radiation 
transport code 0. We considered source cylinders with 
1 m depth in soil, and radii up to 1000 m. The natu¬ 
ral radionuclides were distributed homogeneously within 
the cylinder, and the air and soil compositions and den¬ 
sities followed Eckerman and Ryman m ■ The photon 
emission probabilities and decay branching fractions were 
drawn the recent ICRP publication 107 on nuclear decay 
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FIG. 1. Conversion factors for air kerma per unit concentra¬ 
tion in soil with increasing source radius. 


data [7]- 

The conversion factors increase monotonically in the 
range 1 to 35 m as studied by Gasser et al. pQ. How¬ 
ever, they continue to increase beyond 35 m radius source 
cylinders, and only attain their limiting values to three 
significant figures for source cylinders greater than 700 m 
in radius. The apparent points of inflection on the curves 
in our Fig. |T] are artifacts of the logarithmic abscissa, as 
can be seen from the inset which shows the same data on 
a linear scale. 

Our results yield air kerma conversion coefficients of 
0.042, 0.444 and 0.592 (nGy/h per Bq/kg) for 40 K, the 
238 U series and the 232 Th series, respectively, in the 
asymptotic limit. These values are reasonably consis¬ 
tent with the values of 0.042, 0.463 and 0.604 (nGy/h 
per Bq/kg) respectively found in Saito and Jacob 0j 
(see Table [T] of this letter). Similarly, results recently 
published by Askri m from simulations of large radius 
source cylinders (~ 500 m) were consistent with Saito 
and Jacob |4]. However, the results contrast with those 
of Gasser et al. [I], namely 0.036, 0.357 and 0.482 (nGy/h 
per Bq/kg) respectively for 40 K, the 226 Ra series and the 
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Calculation Source radius Conversion factor (nGy/h per Bq/kg) Nuclear data ref. 


(m) 40 K 238 U series 232 Th series 


Gasser et al. [Tj 

35 

0.036 0.357 a 

0.482 

Nucleide-Lara [6] 

This letter 

35 

0.037 0.383 

0.515 

ICRP 0 

Saito and Jacob [3] 

oo 

0.042 0.463 

0.604 

NEA [H] 

This letter 

(X) 

0.042 0.444 

0.592 

ICRP 0 


TABLE I. Summary of different calculations for air kerma conversion coefficients. 


a Calculation for 226 Ra & progeny only. 


232 Th series. 

For direct comparison we simulated source cylinders 
with radius 35 m. In this case the conversion factors ob¬ 
tained were 0.037, 0.383 and 0.515 (nGy/h per Bq/kg) 
for 40 K, the 238 U series and the 232 Th series, respec¬ 
tively. These values are in reasonable agreement with 
those of Gasser et al. [Tj. The residual differences be¬ 
tween the two estimates for 35 m radius sources, and 
between our latest estimates and Saito and Jacob |3] for 
the true half-space geometry, are of the order of a few 
percent. Therefore most of the ~ 20% reduction in the 
conversion factors obtained by Gasser et al. [1| can be 
explained by a finite size effect. Possible explanations 
for the residual differences once this is accounted for are 
different sources for decay emission data, as suggested by 
Gasser et al. PQ, or other detail variations between the 
calculation methods. These could be differences in the 
energy resolution of the Monte Carlo simulations, sta¬ 


tistical accuracy, or the method used for interpolating 
between the photon fluence rate to air kerma conversion 
coefficients, for instance. 

We encourage Gasser et al. to increase the size of their 
simulations and recalculate the air kerma conversion co¬ 
efficients. Larger radius source cylinders can be simu¬ 
lated quickly and efficiently by transforming the Monte 
Carlo source and detector volumes as outlined by Na- 
mito et al. [12]. Another cheap option is to use a plane 
detector and periodic boundary conditions. Modeling a 
source with a sufficiently large radius will yield results 
closer to the conversion factors established in the liter¬ 
ature. Any remaining differences between their results 
and literature values may then be attributed to the use 
of more up-to-date decay data with lower emission prob¬ 
abilities, and perhaps other methodological differences as 
discussed above. 
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